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FOREWORD 


This report combines and supersedes two previous reports (TPRC Reports 20 
and 22) issued as interim and final reports to an AFOSR grant for studies on substances 
of potential use in the gas laser area. The merging of the two previous reports for 
convenience has necessitated some rearrangement of the text. In addition, a number 
of typographical errors have been corrected. 


The Thermophysical Properties Research Center (TPRC) was founded on 
1 January 1957. Due to the enlargement of the scope of its research operations after 
seventeen years, TPRC and the other affiliated Centers were restructured effective 
1 July 1974 to have an umbrella organization, known as the Center for Information and 
Numerical Data Analysis and Synthesis (CINDAS), to comprise all the Centers, 


The present report, which is Number 29 in the original TPRC Report series, 
is therefore designated CINDAS-TPRC Report 29. 
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ABSTRACT 


Tables of viscosity, thermal conductivity, and diffusion coefficient have been 
prepared for one hundred seven binary mixtures and for nine pure substances, the latter 
being helium, atomic and molecular hydrogen, atomic and molecular deuterium, atomic 
and molecular fluorine, hydrogen fluoride, and deuterium fluoride. Most of the binary 
mixtures are combinations of these substances. The temperature range was from 100 
to 3000 K. For the pure substances, the most recent recommendations were used and 
for the mixtures, the Lennard-Jones 6-12 potential was used for viscosity and diffusion, 
and a linear mixing rule was employed for thermal conductivity. The precision of the 
tabulated values is reviewed, The present state-of-tho-art in predicting these properties 
is briefly surveyed. It is concluded that significant discrepancies can still exist between 
theoretical and experimental quantities. 
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I. INTRODUCTION 


This report, an amalgamation of two earlier reports [1,2], presents calculated 
values of the viscosity, thermal conductivity, and Fickian diffusion coefficient for gases 
and gas mixtures of interest to laser technology. The theoretical background is first 
briefly reviewed in order to acquaint the reader with the basis used in assembling our 
tables and to give him an idea of the assumptions made in the work. This background 
will also enable him to apply our tables in any calculations which he may desire to make. 
The methods which were used to prepare the input to our calculations are then reviewed 
and are followed by the tabulated values. 


Il, THEORETICAL AND EMPIRICAL METHODS 


The three transport properties, viscosity (4), thermal conductivity (k), and 
diffusion coefficient (D) can be defined from the equations 


dm « -pa $8 dt 


which relate the fluxes of momentum, heat and mass to the respective gradients of velocity, 


temperature, and concentration. 


The Chapman-Enskog theory represents a solution of the Boltzmann integro- 
differential equation in which equations for the transport properties result. In CGS [87] 


units, 
jot y = -268.93MT_ (gm om! seo~!) (1) 
o2 Oho (KT /e) 


10? j= AOA, Ty Tht) (cal om-! seo! K~t) (2) 
¢ 


108 p= 2828 FT ochre} (amt neon (3) 


are the three first approximation equations, where T is the absolute temperature and 
the omegas are "collision integrals", Other quantities include M for molecular weight 
and P for pressure. Higher approximations involve additional multiplicative factors 


“ Numbers in square brackets designate references listed in Section V. 


in the above equations, As will be noted below, to the accuracy possible in the present 
work, these higher approximations will be ignored. 


Various limitations exist in the above equations due to assumptions made in their 
derivation, These include (1) only two body collisions ocour, (2) that "classical" rather 
than quantum statistics apply, (3) that all gradients are small, (4) that the theory will 
be applied to systems where the dimensions are much greater than those between particles 
comprising the systems, and (5) that the intermolecular force field is spherically sym- 
metrical, These restrictions, in practice, mean that the theory is usually valid for reason- 
ably dilute gases at normal to high temperatures. The "reasonable" dilution implies 
that our equation of state can be written as PV = RT + By P, while "normal to high" im- 
plies that the temperatures are neither so low nor so high that quantum or dissociation 
effects occur, respectively. ‘The most severe restriction is that the intermolecular force 
field be symmetrical. While, at first thought, this should restrict our attention to the 
monatomic gases, it is found that diatomic and (some) polyatomic gases can also be 
represented In this way. The main modifications which result with polyatomic gases 
are (1) that dipole moments exist (these lead to corrections due to polarity), and (2) 
that internal degrees of freedom contribute to the thermal conductivity, These lead to 
Eucken type corrections due to internal degrees of freedom, (Both corrections are trouble- 
some and will be discussed below.) The basic result theoretically is thus contained in 
the statement that if the law of forces between the molecules in a dilute gas is known, 
the collision integrals (the Q; and Q2. of equations 1-3) can, at least in principle, be 
caloulated and thus their variation with reduced temperature, T* = kT/e (where ¢ is 
the well-depth in the intermoleoular potential function), can be found. Some considerable 
attention was paid in the early stages of this work to reviewing the ability and adequacy 
of present day treatments to evaluate the intermolecular force field function and thus the 
collision integral (omega) functions. Two standard works exist on this subject - the 1964 
"Molecular Theory of Gases and Liquids," reprinted with Notes Added in 1964 [3] and 
the third edition in 1970 of "The Mathematical Theory of Non-Uniform Gases" [4]. As 
might be inferred from the titles, the second citation is more mathematical. It also 
contains a slightly less detailed account of the comparison of theory and experiment. 


Due to the fact that the 1964 updating in the ‘Notes Added’ of [3] is only barely 
adequate - a fact noted elsewhere [5] - the review of theore‘ical methods made here ex- 
tended to a somewhat earlier dats. Essentially, 1954 can be considered as representative. 


The position in 1964 was that the Lennard-Jones 6-12 intermolecular potential 
function had found wide application in calculations of equilibrium and transport properties. 


Ao be nt I a el ll ct 


The indices 6 and 12 refer to the (negative) exponents of the intermolecular separation 
r which have a physical significance in representing attractive and repulsive forces be- 
tween molecules respectively. At that time, reasonably good agreement was thought 

Ms to exist between many of the experimental and theoretical predictions. 


| Since that time, many further studies have been made, In 1963 it was noted (6] 

i that errors existed in the collision integrals of [3] and n retabulation was made. Further 

it recalculations have also been made (see, e.g. [7]) and a later recommendation [8] used 

: the values of [6] at reduced temperatures (T* = kT/e, where ¢ is the well-depth in the 
potential and k is the Boltzmann constant) up to 4.6 and those of [7] at higher temperatures. 
In the present work, our earlier values of collision integrals [6] were used for reduced 
temperatures below 10 and those of [7] above this value. No significant difference exists 

i between the present and [8] choices. The present one was mace solely for ease of 
interpolation, 


a Other molecular models have been used since the publication of the original [3] 
review, An extensive investigation of these was made by Klein and associates [9-11]. 
Whichever model is used, the result is that the representation of both second virial coef- 
ficient and viscosity data over a relatively wide temperature range is equally good for 
many fluids. The range for transport properties is about 2 < T” < 6 where the ¢/k in 
the T" refers to the value determined using the 6-12 parameter. The adjective "about" 
refers to the fact that available viscosity data are imprecise - a statement that will be 
amplified during the discussion of the tabular values below. 


The temperature range of interest to the present program involves reduced 
temperatures (on the Lennard-Jones 6-12 scale) of 0.4 to 70, depending on substance, 
Particularly for the light gases and mixtures, reduced temperatures above 6 occur, 
In this region the energies of intermolecular collisions may be large enough that the 
attractive forces play a relatively weak role. This is confirmed by an examination of 
the tables of Klein [7] for the Lennard-Jones m-6 potential, where, in addition to the i 
’ value of m = 12 previoualy discussed, Klein tabulates collision integrals Q4, and Qo. 
4 (and other functions) for m #9, 15, 18, 21, 30, 60, and 75, If one plots the collision 
- integral Qh. as a function of T* on logarithmic paper, one can conclude that the repre~ 
sentation Qo) = a(T) "isa quite acourate approximation for reduced temperatures above 
about 15. Furthermore, one can conclude that a definite relationship exists between ’ 
the nandm, Finally, that the rigid sphere repulsion force model is approached as m : 
{noreases, These resulta check with the assumption [3,4] that y= AT® if the repulsive 4 
potential only exists. 
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The m-6 model was then refined by Hanley and Klein [12,13] to the m-6-8 model, 
where the 8 refers to the dispersion forces. (For a description of these, see, e.g. page 
29 of [3].) The initial expectation was that the coefficients A; and A, in the equation 


E(r) sAgr?+ Agrt+ BL ™ (4) 


in this m-6-8 potential should be reconcilable with theoretical estimates, However, 
the results [13] of comparing theory with experiment have proved alesepointios: The 
following extracts from [13] may prove informative. 


1 (page 1747): Any spherical potential oan never properly depict the interaction 
between two polyatomic molecules - as the molecules collide a number of 
orientation-dependent force laws are followed - one for each relative orientation. 


2 (page 1747): Collisions between polyatomic molecules can be inelastic, The 
effect of the molecular internal degrees of freedom has then to be considered. 


3: In addition, the discussion on pages 1750 and 1751 shows that the coefficients 
A, and A, deduced from theory and experiment agree poorly. On page 1761 
it ia stated that "adding more parameters to our potential for additional flex- 
ibility would not necessarily reduce the disagreement with theory even when 
additional parameter are added to the most uncertain part of the potential, 
the repulsive branch." 


Difficulties also exist for fluids for which several transport properties have 
supposedly been accurately measured. Hanley, McCarty, and Intemann [14] were unable 
to reconoile extremely carefully calculated and measured thermal conductivities and 
viscosities for hydrogen. This was also noted for inert gasea by Kestin (15]. Both workers 
favored the viscosity data, supposing an (uncertain) error to exist in the thermal con- 
ductivities. However, both workers relied on the basic accuracy of the Chapman-Enskog 
theory, Another difficulty was noted by Dahler [16] who said that there was a remote 
possibility of failure due to this theory. A stronger viewpoint is taken by Brush [17] 
who decided that the interatomic force concept was inadequate. This is a conclusion 
backed by at least plausible arguments which he makes - one of which ia identical to the 
first cited remarks of [13] above. 


All of thege results are disturbing to a short term program to caloulate transport 
properties, In addition to the above the writer consulted [18-22] and many other leas 
detailed reviews in whioh further studies of interatomic and intermolecular potentials 
result, Another viewpoint can be taken which is soientifically more pessimistic but which 


can yield results of somewhat predictable accuracy. That is, choose a potential for which 


the collision integrals and related functions are reasonably accurate and available. Fur- 


thermore, use reasonably accurate parameters in the sense of Klein (8]. This should 
seoure at least reasonable results in the reduced temperature range 2 < T <5 and pos- 
sibly for 2 <T* <10. While the accuracy at higher temperatures may be marginal the 
hope ia to adjust the fit by varying the repulsive index. 


In this work, following the survey which led to the conclusions above, the Lennard- 
Jones 6-12 potential was selected for generation of the tabular values. While some m-6 
fittings for simple gases have been made, the availability of data could not in the time 
available lead to acourate values of m for mixtures, The values of the collision integrals - 
and the fittings - for the m-6-8 potential have not yet been published", However, some 
predictions of the third virial coefficient using it were in disappointingly poor agreement 
with experiment [24] . 


One oan conclude that no intermolecular potential is oapable of representing the 
data for all properties to within the experimental error. Thia is due to orientation effects 
and inelasticity in collisions coupled with the inability to perform the mathematical cal- 
culations to sufficient accuracy and the fact that high energy collision and spectroscopic 
data may yield still different intermolecular potential values. This conclusion is in part 
supported by a subsoquent analysis [54,65] of the ability of the m-6-8 potential to simul- 
taneously fit viscosit,; and second virial coefficient data for non-polar polyatomic gases, 
It was found that additional terms were necessary to represent the virial coefficient be-~ 
havior, This implied that a more comprehensive potential should then likewise be used 
to caloulate the viscosity (and, of course, other properties). This is an enormous com- 
plication, both from the fundamental diffioulty of calculating functions for comprehensive 
potentials and from the more practical viewpoint of fitting such potentials to the exper- 
imental data, While it is atated [66] that "work is in progress to modify the m-6-8 
potential to include non-spherical effects" one reads on the same page that "one assumes 
that the molecules collide with a fixed relative orientation and collision integrals are 
caloulated for each orientation independently, The final values of the collision integrals 
are then determined by a statistical averaging over all these possible orientations". 
Here is the important fact that the caloulations to date assume a fixed relative orienta- 
tion. This, coupled with the averaging requirement will, in the opinion of the writer, 
introduce sufficient uncertainty so as to negate possible improvements by modifying the 
potential. Indeed, it appears probable that this stage has already been reached from 
The report "Tables of Collision integrals and Second Virial Coefficients for the (m, 6, 8) 


Intermolecular Potential Function" by Klein, Hanley, Smith, and Holland has sinoe 
(June 1974) been published as NSRDS-NBS-47, 157 pp., 1974. 
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a fundamental viewpoint by the failure of the potential to correctly model the dispersion 
coefficients, as already mentioned in the discussion of Eq. (4). An extensive review 

(22] of some aspects of intermolecular forces does nothing to dispel the conclusion that 
realistic potentials are not yet available for molecules of complex structure. As remarked 
by Svehla [23], unless a consistent potential is selected, values of combinations (such 

as the Prandtl Number) of properties may vary irregularly from substance to substance, 
This is even more important in preparing tables for mixtures for which the composition 
varies. Hence, in the present work, the Lennard-Jones 6-12 potential was selected as 

the general basis for generating the tabular values rather than a more sophisticated model, 
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a Ill, EQUATIONS FOR CALCULATING TRANSPORT PROPERTIES OF 
4 PURE GASES AND GAS MIXTURES 


1. GENERAL DISCUSSION 


As noted above, Eqs. (1-3) represent first approximation equations for the 
variation of these properties with temperature for pure gases. In order to use them, 
the collision diameter, o, must be known. In addition, the well depth ¢/k is needed so 
that the collision integral Qi; or {Q_ can be determined for the particular reduced tem- 
perature IT of kT/e of interest. 


The result of detailed caloulations has been to reproduce tables of N,, and/or 
Qgq as a function of T*, An early example appears in [3]. Later tables were intercom- 
pared in (6] and again in [7]. The reader is referred to [6] for the extensive (3,320 
values) tables of ,; and Qo. used in this work for reduced temperatures below 10. Table 
117 contains original and interpolated values of [7] for Qi, Qoo, and A, all as functions 
of kT/¢ for the Lennard-Jones 6-12 potential, (The A function is described in the sec- 
tion on viscosity of mixtures.) While various attempts have been made to represent the 
variation of these functions with reduced temperature analytically, most of these have 
suffered from the difficulty that the tables of [3] were used ag inputs to the analyses, | 
Ag noted elsewhere [6-8], the [3] tables do contain relatively large errors. Thus, the { 
analyses need repeating using more accurate input. While two recent studies have ap- 
peared [8,25] which use the [7] tables as input, it was not considered feasible to program 
these in time for the present report, The more recent fits [8,25] also are limited to 
the (12,6) potential. It is considered desirable to derive a convenient representation i 
for the (m,6) or, even better, for the (m-6-8) model. To allow for the polar-nonpolar 

4 


interactions, possibly a (3-m-6-8) function should be devised. 
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2. VISCOSITY 


Lacking such input at the present {ime, the omega tables of [6] and [7] were used 
numerically, in the ranges as described above, with Eqs. (1-3)", to generate the present 
tables. The parameters o and ¢/k were taken from the work of Svehla [23], This is 
the main previous effort to present a unified body of thermal conductivity and viscosity 
data based on a common potential and, furthermore, which selects molecular parameters 
which are reasonably self-consistent, By this is meant the fact that the estimation methods 
for compounds result in parameters which yield reasonable agreement with experimental 
: data for compounds. The adjective reasonable is used to note that better agreement for 
2 specific properties and substances has been obtained by using different potentials. Pre- 
sently lacking anywhere in the literature is a complete across-the-board reanalysis, 
except possibly for the inert gases [e.g. 16], The comparisons for thermal conductivity 
in [26-28] do reveal that the Svehla values may be in error by possibly up to twice the 
accepted limits for some substances, For other materials, a very reasonable fit occurs 
in general, 
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For interpolation purposes, linear interpolation of the prasent tables is considered 
perfectly acceptable. If an empirical representation is needed, it is suggested that the 
Q functions of Eqs. (1-3) be replaced by 


Q=a+b/T +0/T? +d/T3 (5) 


Aa will be evident from [28-30], use of Fq. (6) to represent viscosity and thermal con- 
ductivity results in negligible errors over large temperature ranges, even using d = 0 
in many cases, 


For mixtures, the equations are considerably more complex. For gaseous viscosity 
the mixture viscosity, bay is given by 


142 
MM = XY 8) 
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In these equations M,;, My, x,;, and x, are the molecular weights and mole fractions of 
components 1 and 2 and the wij, Wie, ANd pop represent the viscosities of component 1, 
a hypothetical gas 12 and component 2, The gas 12 is assumed to obey the same poten- 
tial as gases 1 and 2 but to have a molecular weight of 2M, M,/(M, + Mg), @ collision 
diameter of (0; +0,)/2 and a well depth ¢/k of (¢;/k) (e9/k). Aj, ig @ collision inte- 
gral function [3], tabulated simply as A in Table 117. 


In the various tables, the three columns under viscosity are values of wii, jie» 
and jigq for the different temperatures listed. As shown by Table 117, A varies from 
1.1106 at a reduoed temperature of 10 to 1.1323 at T*# = 100, or 1.1214 + 0,0109 over 
this temperature range. To a better precision, the approximation 
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A = 1,09340 +0.020 log T# (7) 


will represent the [7] A values for 15 < T* < 100 to an average error of 2x 10" anda 
maximum error of 1x 1073 in A. This approximation may be found useful in program- 
ming Eq. (6). 
For multicomponent mixtures and for higher approximations, the careful 
comparisons of theory and experiment by Brokaw, et al. [31,32] may be found useful. 
The lack of accurate experimental data for most of the mixtures considered here was 
| one reason why the higher approximation equations were not considered further in the 
" present work, 


38. THERMAL CONDUCTIVITY 


Equation (2) above for thermal conductivity ig only found to be applicable to the 
caloulation of thermal conductivity for monatomic species, where no contributions of 
internal degrees of freedom to the thermal conductivities ocour. The classical paper 
on the subject since [3] is that of Mason and Monchick [33] where polyatomio and polar 
gases were considered, In addition, Chapter 13 of [4] provides an excellent review, 
as do publications by Brokaw [34] and by Touloukian, Liley, and Saxena [27]. The equa- 
tions involved are complex and the reader is referred to the cited sources for details. 
The conclusion one can draw from a study of these works is that even an extremely thor- 
ough and painstaking analysis can result in thermal conductivity values and/or the related 
Eucken factors that disagree significantly with experiment. The publications by Hanley, 
et al. [14] and Kestin [15] may also be olted here, Further evidence of disagreement 
may be found in [33,38]. For many substances, simple equations representing the ther- 
mal conductivity aa linear and quadratic functions of temperature have been found adequate 
[28-30, 36,37] aa has also the representation of Eq. (2) where the omega function 
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approximation of Eq, (5) is used [28-30]. The conclusion to be drawn from this is that 
the relative magnitude of the internal contributions to the total energy decreases to about 
a constant value at high temperatures, While it cannot be neglected, the approximation 
of Svehla [23] oan be seen from the comparisons in [28-30] to yield reasonably accurate 
results at high temperatures, though not at low temperatures. Hence, in this work, 
values based essentially on experimental data were used for the elements and simple 
compounds and were taken from [26-29]. For other compounds the values of [23] were 
used, or analogous estimation methoda, 


The computer print-out was, in fact, programmed so that many of the tables of 


‘yiseosity and diffusion coefficient appearing in this report are edited copies of the com- 


puter sheets, The thermal conductivity values were based on experimental data, cor- 
related with theory if possible as described in [29], or were generated by applying 
correctiona for internal degrees of freedom to the calculated translational effect, In 

the latter case, some substances showed a significant divergence between the theoretical 
and the experimental values. In such cases, the caloulated values were reduced to bring 
them into general accord with the experimental data. In our programming we adopted the 
fittings of Neufeld, Janzen, and Aziz [25] to generate the 6-14 collision integrals. It 
should be noted that their representation is only valid for reduced temperatures from 0.3 
to 100 and not from 0.1 to 0.3 or from 100 to 400. Some problema were encountered be- 
fore this fact was realized. While the lower end gives no practioal difficulty, substances 
with low characteristic well depths, ¢/k, oan yield reduced temperatures over 100. A 
suitable blocking procedure was thus written into the program to prevent the computer 
extrapolating the calculation into regions where the reduced temperature exceeded one 
hundred. For the few cases where this occurred, tables were numerioally generated 
using the approximations 


Oh, = 1.10670 (TH) “0 1675 
(4) 

Qhy = 0.07188 (TH) “He 1825 
whioh have been found to yield collision integrals acourate to 1/4 percent for reduced 
temperatures over 28, In fact, the use of the Lennard-Jones potential at such high re- 
duced temperatures is physically unrealistio due to the predominance of the repulsive 
part of the potential for which an exponential dependence on intermolecular separation 
is preferable to a power one, although, as mentioned earlier, neither representation 


ja exact. 


For mixtures, the equations are even more complicated and the analysis of [27] 
centered on approximations for the functions equivalent to the Aj, of Eq. (6). A careful 
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study was made of the thermal conductivity - composition curves of [27] but no simple 
trends were evident. A cogent comment was made by Lehmann [38] and cited by 
Bretsznajder [89] that simple linear additivity is as reasonable as more sophisticated 
caloulations for binary mixtures for which the molecular weight of the pure components 
do not differ by more than a factor of 4. Examination of the graphs of [27] reveal that - 
in a surpriaingly large number of cases ~ the variation with composition is easentially 
linear, even for some mixtures where exceptions might be expected to ocour. A more 
accurate caloulation would involve knowledge of relaxation times, some of which are 
now available under an AFOSR contract [40], In any future work on these substances, 
or those of lower priority, it is considered highly desirable to attempt to establiah at 
least a function correlation of Eucken factor with molecular weight and/or chemical 
atructure, 

The thermal conductivity columns in the tables thus present values for the pure 
component and for the equimolar mixture, tn :\ny future work it is only antioipated that 
the equimolar values may show significant chunges. 


4. DIFFUSION COEFFICIENT 

For diffusion coefficient values, Eq. (3) may be used for the self-diffusion 
coefficient and for interdiffusion, the equation is similar to (3) except that the o and 
the ¢ now refer to a hypothetical gas 12, just as for viscouwity. No attempt was made 
to introduce a variable coefficient into Eq. (3) as suggested by Wilke and Lee [41], 


The three diffusion coefficient columns, tabulated for atmospheric pressure, list 
self-diffusivity for the pure components and the interdiffusion for the binary pair, 


IV, REVIEW OF INFORMATION AVAILABLE ON TRANSPORT PROPERTIES 
FOR SUBSTANCES OF INTEREST AND OF METHONS USED TO PREPARE 
TABLES OF RECOMMENDED VALUES 


1, DEUTERIUM 


The Lennard-Jones potential and parameters of Svehla [23] were used to caloulate 
the viscosity and diffusion coefficient values. The thermal conductivity values to 400 K 
[20,83] were extrapolated using a log-log plot. 


2. HELIUM 


Until recently, up-to-date information on the properties of hellum has boen 
nonexistent. Within the last few years the position has ohanged for the better, 
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For the viscosity, the tables of Maitland and Smith [42], Watson [43], and 
Petersen [44] are in very good agreement, In somewhat less good agreement are the 
Tsederberg, et al, tables [45] while Angerhofer and Hanley [46] present tables for rel- 
atively low temperatures, The values presented here are essentlally a mean of the [42-44] 
sets for temperatures to 2000 K, Values at higher temperatures were derived from an 
analysis of the lower temperature values which showed that p ™ aT™ where n> 0,702, 

The tabular values Were actually derived by assuming n to be a weakly varying linear 
function of temperaiure. The agreement between our tabulated values and any of the 
source values listed above is well within five percent. The absolute accuracy may, in 
fact, be a little better than thie figure, 
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Thermal conductivity values to 2500 K have been tabulated by the author elsewhore 
[28,29]. The present tables are merely a reproduction of these with a finite difference 
extrapolation to 3000 K. Irom 100 to 700 K the acouracy should be two percent, from 
700 to 1500 K five percent, and above 1500 K ten percent. 


Values of the diffusion coefficient were obtained using the equation suggested in 
the Mason and Marrero [47] compilation for ‘Ho-‘He diffusion. Values tabulated in 
Landolt-Bornstein [48] for helium agree to within about +10 percent. 


8. ATOMIC HYDROGEN 


To ensure better consistency in the mixture tables, the tables of viscosity and 
thermal conductivity of atomio hydrogen tabulated by Svehla [23] were used even though 
some disagreement with the later Hanley, et al. [14] values ocourred. The disagree- 
ment amounts to less than six percent In viscosity and eight percent in thermal conduc- 
tivity. The diffusion coefficients, obtained from the work of Mata [62] were extrapolated 
above 2000 K assuming a logarithmic dependence on temperature. Possibly a ten percent 
uncertainty is reliable. 
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4. HYDROGEN 


Three seta of tables for viscosity, by Maitland and Smith [42], Watson [43], and 
Hanley, McCarty, and Intemann [14] show insignificant differences between 100 and 2000 K. 
These were extrapolated to 3000 K to parallel the lower earlier tables of Bambanek [48] 
and Svehla [23]. An error of a few peroent below 1000 K increasing to 6-10 percent at 
3000 K appears possible, 


The thermal conductivity values recommenced in [28,29] to 1000 K were extrapolated 
to 8000 K to merge with the earlier values of [23]. The uncertainty should be a few 
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percent below 500 K, this increasing to possibly ten percent at the highest temperatures 
tabulated, 

A plot of the self diffusion coefficient values tabulated by D'Ans-Lax [49] and 
in Landolt-Bornstein [50] revealed an almost linear variation of log D with log T. The 
tabulated values, obtained from this graph (extrapolated where necessary), should be 
accurate to five percent below 500 K, the uncertainty increasing to ten to fifteen percent 
at 3000 K. 
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§, OTHER PURE SUBSTANCES 


While the properties of helium, hydrogen, and deuterium discussed above were 
generated by the various methods there described, for other pure substances such as 
argon, fluorine, and nitrogen the parameters of Svehla [23] were used with the Lennard- 
Jones 6-12 potential firat to generate the tabular values in order %o achieve the maximum 
consiatency in the tables for the mixtures. Parameters for substances not considered 
by Svehla were estimated uaing procedures identical or similar to those that he used, 

A collection of the parameters obtained in this way appears in Table 118. 


No table of intermolecular parameters for all the gas pairs considered here is 
given as many ofthetables* were generated from a computer program in whioh the mix- 
ture parameters were internally generated from the pure component values. The usual 
[3,86] methods for generating the parameters for mixtures were followed. 
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6. MIXTURES 


For all mixtures, the viscosities tabulated under the mixed heading are the 
viscosities of a hypothetioal mixed gas, as desoribed in the text following Eg. (6) where 
the parametera used in the potential are obtained using the appropriate combination rules 
with the pure substance parameters of Table 118. The thermal conductivities tabulated 
under the mixed heading are for an equimolar mixture and assume a linear dependence 
on composition while the diffusion coefficients are for a hypothetical mixed gas. 
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Experimental data are laoking for many of the systems of mixtures, The error 
in our tables can be assumed, in general, to be about twenty percent greater than the 
larger of the two uncertainties for the pure component viscosities or diffusion coefficients, : 
For thermal oonduotivities for systems for which the pure component molocular weights 
differ by lesa than a factor of four this should also apply. For thermal conduotivity of i 
systema of more disparate molecular weights, greater error may ocour. 


¥ Minor differences exist ina (ew cases between values of properties for the same substanco 
tabulated in [1] and [2] and thus also in the different tables reproduced here, 
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minded that appropriate conversion factors are needed to convert from CQS to SI 
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in volumes 11 and 3 of the TPRC Data Series. 
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TABLE 117. COLLISION INTEGRAL AND ANOTHER FUNCTION FOR THE 
LENNARD-JONES 6-12 POTENTIAL FUNCTION 


On Che A T™ hy Ch A 
0.7422 0.8242 1,1105 55 0.5674 0.6405 1.1287 
0.7806 0.8124 1.1119 56 0.6658 0.6387 1,1288 
0.7202 0.8017 1.1131 57 0.5642 0.6870 1.1289 
0.7109 0.7921 1.1142 58 0.5627 0.6853 1.1291 
0.7026 0.7884 1.1152 69 0.5612 0.6837 1.1292 
0.6948 0.77684 1.1161 60 0.5598 0.6322 1.1293 
0.6877 0.7681 1.1169 61 0.6584 0.6806 1.1294 
0.6811 0.7612 1.1176 62 0.5570 0.6290 1.1295 
0.6750 0.7549 1.1188 63 0.6566 0.6275 1.1297 
0.6698 0.7489 1.1190 64 0.6642 0.6260 1.1298 
0.6640 0,7483 1.1198 65 0.5528 0.6246 1.1209 
0.6698 0.7383 1.1201 66 0.5514 0.6282 1.1800 
0.6642 0.7383 1.1206 67 0.5501 0.6218 1,1301 
0.6496 0.7288 1.1211 68 0.5489 0.6205 1.1301 
0.6453 0.7297 1.1216 69 0.56477 0.6191 1.1802 
0.6413 0.7194 1.1219 70 0.6465 0.6177 1,1803 
0.6974 0.7168 1.1223 71 0.6463 0.6164 1,19304 
0.6886 0.7118 1.1227 72 0.6441 0.6151 1,1906 
0.6301 0.7076 1.12381 73 0.6429 0.6189 1.1305 
0.6267 0.7040 1.1234 74 0.5417 0,6126 1.1306 
0.62384 0.7008 1.1237 76 0.6406 0.6118 1,1307 
0.6203 0.6972 1.1240 76 0.6306 0.6101 1.1308 
0.6172 0,6930 1.1244 77 0.5384 00,6089 1.1309 
0.61438 0.6908 1.1247 78 0.6873 0.6077 1.1909 
0,6114 0.6878 1.1249 79 0.6362 0.6065 1.1310 
0.6086 0.6849 1.1252 80 0.6962 0.6053 1.1311 
0.6060 0.6820 1.1264 81 0,6341 0.6042 1,19312 
0.6036 0.6793 1.1267 82 0.6391 0.6031 1.1913 
0.6010 0,6766 1.1269 83 0.53821 0.6020 1,1813 
0.6986 0.6741 1,1261 84 0.53811 0.6009 1.1814 
0.6962 0.6715 1.1263 85 0.5301 0.5998 1.1915 
0.5840 0.6690 1.1265 86 0.5291 0.6988 1.1316 
0.5919 0.6665 1.1267 87 0.6282 0.6077 1,1316 
0.5897 0.6640 1.1268 88 0.56272 0.69614 1,1317 
0.5878 0.6620 1.1270 89 0.6263 0.6056 1.1317 
0.5854 0.6599 1.1272 90 0.5254 0.65946 1.1318 
0.6835 0.6578 1.1274 91 0.5245 0.6936 1.1318 
0.5816 0.6557 1.1276 92 0.5236 0.5926 1,1318 
0.5797 0.6587 1.1277 93 0.6227 0.5017 1.1319 
0.6778 0.6517 1.1279 04 0.6218 0.6907 1.1319 | 
0.5758 0.6490 1.1281 95 0.5209 0.6897 1.1920 
0.5741 0.6477 1.1282 96 0.5200 0.5888 1.1920 
0.57239 0.6468 1.1283 97 0.6191 0.5879 1.1921 
0.6706 0.6440 1.1285 98 0.5183 0.6869 1.1921 
0.5690 0.6422 1.1286 99 0.5175 0.5860 1.1822 


| 
100 (0.5167) 0.58611. 1828 | 
| 
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TABLE 118. LENNARD-JONES 6-12 PARAMETERS FOR PURE GASES 


Gas e/K(°K) ~—g (A) 1//M 
He 10,22 2.561 0.48980 
Ny 113. 50 3.566 0, 18894 
A 138,20 3, 287 0, 16822 
Hy 59.70 2,827 0.70430 
D 36.20 2.664 0.70712 
NH, 558. 30 2.800 0. 24232 
ND, 3387.00 2,900 0. 22367 
NgH, 760.00 3,180 0.17666 
NoD, 462.00 3. 380 0, 16664 
Fy 112.60 3, 307 0.16223 
F 112.60 2,968 0, 22942 
NI, 175.00 4,154 0.11868 
NF, 240,00 3, 880 0, 08805 
CIF, 448,30 4.900 0.08770 
HF 330.00 3.148 0.22367 
DF 198.10 2, 826 0.21223 
HCl 344.70 3, 939 0.16658 
DCl 208.00 2,980 0.16340 
H 37. 0U 2,701 0.99210 


D, 35.21 2.982 0, 60000 


